Finite Conductivity Minimum in Bilayer Graphene without Charge Inhomogeneities 
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Boltzmann transport theory fails near the linear band-crossing of single-layer graphene and near 
the quadratic band-crossing of bilayer graphene. We report on a numerical study which assesses the 
role of inter-band coherence in transport when the Fermi level lies near the band-crossing energy 
of bilayer graphene. We find that interband coherence enhances conduction, and that it plays an 
essential role in bilayer graphene's minimum conductivity phenomena. This behavior is qualita- 
tively captured by an approximate theory which treats inter-band coherence in a relaxation-time 
approximation. On the basis of this short-range-disorder model study, we conclude that electron- 
hole puddle formation is not a necessary condition for finite conductivity in bilayer graphene at zero 
average carrier density. 



INTRODUCTION 



The robust conductivity of nearly neutral graphene 
sheets 1 -^ is interesting from a theoretical point of view, 
awkward 3 for some potential applications, and among the 
most unexpected of graphene transport study discoveries. 
As a function of ambipolar carrier density the minimum 
conductivity is ~ e 2 //i, with relatively small sample to 
sample variation. The generally accepted explanation^— 
for this property starts by recognizing the influence of 
randomly distributed charged- impuritie s 13 ' 14 which in- 
duce electron-hole puddle s 15 ' 16 in graphene when the 
global average carrier density is low. Partly because oS 5 - 
the role of Klein tunneling in Dirac-like systems, a net- 
work of conducting puddles can account for global con- 
duction that remains finite when the average carrier con- 
centration falls to zero. There are, however, indications 
that this explanation is incomplete. In particular, sus- 
pended graphene 1 ^— samples still exhibit a minimum 
conductivity even though charged impurities appear to 
be removed upon annealing and puddle formation should 
therefore be suppressed. The present work is motivated 
by the view that graphene's minimum conductivity phe- 
nomena are more general than sometimes thought, and 
not necessarily associated with smooth inhomogeneities. 

Independent of disorder character, transport near the 
band-crossing energies of graphene systems differs from 
transport near typical semiconductor band extrema in 
three important ways: i) the absence of an energy gap 
between the conduction and valence bands, ii) the pecu- 
liar momentum-dependence of inter-sublattice hopping in 
graphene systems that leads to the Dirac-like electronic 
structure and iii) in the case of single-layer graphene the 
linear band dispersion which causes the two-dimensional 
density-of-states to vanish in the absence of disorder. The 
goal of this paper is to shed light on which of these as- 
pects is responsible for conductivity minimum phenom- 
ena. Since experiment indicates that there is no essential 
difference between the minimum conductivity behavior of 
single and bilayer cases, the dispersion law does not ap- 
pear to play an essential role. The minimum conductiv- 



ity is also finite in suspended bilayer graphene^ samples, 
even though the charge carriers in this case exhibit the 
same parabolic 2 — dispersion that is found in conventional 
two-dimensional electron systems. We therefore focus on 
bilayers, and on the role of momcntum-sublattice cou- 
pling in the absence of an energy gap. This problem has 
received relatively little theoretical attention 2 ^—. 

Momentum-sublattice coupling in bilayers is well de- 
scribed by the 7r-band envelope-function effective band 
Hamiltoniar>ii2i 
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Here m ~ 0.05m is the effective mass, mo is the bare 
electron mass, k is the two-component particle momen- 
tum, and the matrix structure originates from the layer 
and sublattice degrees of freedom. The Hamiltonian Hq 
does not contain the trigonal warping term and just rep- 
resents the minimal model where the conductivity min- 
imum does not vanish. The spectrum of Hq consists of 
parabolic conduction and valence bands that touch at 
eigenenergy E = 0. The sublattice degree-of- freedom 
is frequently viewed as a pseudospin in order to exploit 
analogies between spin-orbit and pseudospin-orbit cou- 
pling. From this point of view Hq can be considered as 
expressing an effective Zeeman coupling to pseudospins 
that has a strength hflk = h 2 k 2 /m which is momentum- 
magnitude dependent, and a x-y plane orientation an- 
gle <f> — 20k where 0k is the two-dimensional momen- 
tum direction. The pseudospin precession axis therefore 
changes whenever an electron is scattered between mo- 
mentum states. When the precession frequency VLk is 
larger than the momentum scattering rate r _1 , the pseu- 
dospin precesses a few times between collisions and any 
initial transverse component is likely to be randomized. 
The conductance minima phenomena occur for energies 
E near zero for which f^r is always small and pseu- 
dospin components transverse to the precession axis are 
not expected to randomize. This observation alone sug- 
gests the possibility that atypical quantum effects could 
play a role. This is what we can see in Fig. [TJ The con- 
ductivity never falls to zero for any reasonable choice of 
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FIG. 1: The dotted curves depict the electrical conductivity 
of bilayer graphene (per spin/valley) as a function of carrier 
concentration computed according to the Kubo formula (|2| 
for the series of model parameters specified in Table [T] The 
solid lines correspond to the analytical approximation which 
is the sum of the Drude conductivity ctd and an interband 
coherent correction Act given by Eq. J5|. The inset illustrates 
the decomposition of the conductivity for disorder model A 
into intra- and interband coherent contributions proportional 
respectively to the intra- and interband terms in the velocity 
operator in Eq ([2]). 

parameters as long as the interband coherence is included 
in the model, even when charged impurities are absent 
and electron-hole puddle formation^^ is not expected. 
We focus solely on the zero temperature limit. The fi- 
nite temperature ^ 1 1 1 3 1 can lead to the thermally excited 
carriers which may spoil the interband coherence effect. 
The intervalley scattering is also assumed to be absent 
here. 



II. KUBO AND BOLTZMANN THEORIES 

We have evaluated the conductivity numerically us- 
ing the non-interacting particle Kubo formula expression. 
This approach has the advantage that it is exacts, or at 
least would be if computational resources were infinite. 
On the other hand it does not lend itself to a satisfying 
qualitative understanding. We therefore compare our nu- 
merical results with those predicted by a heuristic semi- 
classical theor y 33 ' 34 that captures inter-band coherence 
corrections to the Boltzmann equation. We first com- 
ment briefly on these two approaches. 

The finite-size Kubo formula for the static conductivity 

is, 

ihe 2 \ ^ Jg n - ! { E n , (n\v x \n')(n'\v x \n) 
aK ~ 1? ^ K - E n > E n -E n ,+i V ' [ ' 

n,n' 

where v is the velocity operator, is the Fermi- 

Dirac distribution function, and \n) denotes an exact 
eigenstate of the Schrodinger equation for a finite-size 
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0.30 
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0.81 


0.13 
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0.25 


0.83 


0.97 


0.10 


C 


0.20 


0.66 


1.22 


0.08 


D 


0.15 


0.50 


1.62 


0.06 



TABLE I: Parameters for Fig. [TJ r is the momentum relax- 
ation time, fi — er/m is the mobility of carriers, n s is the 
concentration of short range scatterers with strength fixed at 
a value uo = it 2 ft 2 /5m small enough to validate the golden- 
rule life-time expression, and 8E — 2nh 2 / L 2 m is the level 
spacing for sample size L = 1.8 x 10 -5 cm. At this sample 
size dependence on L is weak. The momentum cut-off k* and 
L fix the Hamiltonian matrix dimension at 3362 x 3362. The 
computations have been performed at zero temperature. 



disordered system with periodic boundary conditions: 
(H + U)i> n = E n ^ n with U(r) = u £f " 5(r - R*) for 
the short-range disorder model we consider. The scat- 
tering locations and potential signs are random. We 
solve the Schrodinger equation using a large momentum- 
space cutoff k* w \/5 • 10 13 cm -1 which corresponds to 
the energy scale at which the split-off bands of bilayer 
graphene become relevant and our two-band model no- 
longer applies. 

The physical dc conductivity can be obtained from 
Eq. ^ by extracting the limit in which the system size 
first approaches oo, and then r\ approaches zero maintain- 
ing a value larger than the typical level spacing 8E. For 
the model considered here 5E — (2nh 2 ) / '(mL 2 ) where L 2 
is the finite-size system area. The finite value of rj can be 
understood as representing energy uncertainty due to the 
finite lifetime of electrons in a system coupled to source 
and drain reservoirs. To eliminate the influence of the 
bath on the conductivity itself, the momentum relaxation 
time r due to internal scatterers must be much smaller 
than h/rf&. We estimate r using the Fermi golden-rule 
expression: r = 2h 3 jran s u\ where n s = N s /L 2 is the im- 
purity density. Since the smallest possible SE is limited 
by numerical practicalities, we can estimate the conduc- 
tivity only for relatively strong disorder. Conductivities 
obtained directly from Eq. @ undergo the phase coher- 
ent fluctuations; we simulate macroscopic system con- 
ductivities by averaging the conductivity over an energy 
interval containing 10-100 levels, over boundary condi- 
tions, and over several disorder potential realizations. 
Note, that the conductivity fluctuation amplitude turns 
out to be essentially smaller that e 2 /h near the neutral- 
ity point. This makes our numerical approach reliable 
for the conductivity minimum evaluation. 

Below we compare our numerical results for the con- 
ductivity to an analytic modified Boltzmann equation 
theory. When coherence effects are retained the distri- 
bution function /(k) becomes a 2 x 2 matrix with band 
labels. The steady state limit of its equation of motion 
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i{H J(k)} =I[/(k)], 



(3) 



where E is an electric field small enough to justify linear- 
response theory, I[/(k)] is the collision integral which 
accounts for disorder scattering, and the commutator 
[H ,f(k)} accounts for the difference in time evolution 
between conduction and valence band eigenstates. When 
the collision integral is evaluated to leading (second) or- 
der in the (configuration averaged) impurity potential, 
the collision term (including its off-shell term o 34 ^ 36 ) re- 
duces to the simple matrix relation-time form, I[/(k)] — > 
— /W(k)/r, where /W is the deviation from equilibrium. 
This is a remarkable property of the two-band bilayer 
model with <5-function scatterers. In the Ho eigenstate 
basis, the density-matrix linear response then reads 
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Here, Ek± are the eigenvalues of Hq, f2fc = hk 2 /m, and 
v Jjff / is the velocity operator written in the Hq eigenstate 
basis. Given this approximation for the linear response 
of the distribution function, it is easy to calculate the 
Boltzmann conductivity: ctb = jx/E x where j x is the 
electrical current, j = ej Air [v/W(k)]. Note that 

neither v nor / W(k) are diagonal, and that j therefore in- 
cludes interband coherence contributions. The intraband 
contribution to the conductivity stems from the diagonal 
terms in Eq. (j4]) and is given by the simple Drude for- 
mula <td = e 2 nr/m, where n is the carrier concentration 
n = kp/(4ir) with kp being the Fermi momentum. 



III. RESULTS 
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FIG. 2: Comparison between Kubo conductivities ([2]) of (a) 
the decoupled band model and (b) bilayer graphene. These 
results were obtained for a series of models with identical 
golden-rule relaxation times r = 0.3 ■ 10~ 13 s, and sample 
sizes L — 1.8 x 10 -5 cm. (The concentration of scatterers n s 
was adjusted appropriately in each case.) One can see that the 
conductivity minimum for the decoupled band model vanishes 
whereas for the bilayer model it is finite and insensitive to the 
scattering potential strength. The thick solid lines show the 
naive prediction of (a) Drude theory and (b) our interband 
coherent Boltzmann model with golden-rule relaxation times. 



Numerical results for the dependence of Kubo conduc- 
tivity on carrier density are presented in Fig.[T]for a series 
of model parameter values summarized in Table F Our 
main finding is that the conductivity remains finite as 
the carrier density approaches zero. We do not observe 
any systematic dependence of the minimum conductiv- 
ity, CT m in ~ 0.7e 2 /h per spin and valley, on model system 
parameters. 

There are two elements in our model which couple the 
two bands in the Hamiltonian ([1} and both are important 
for the conductivity minimum phenomena. The first is 
the velocity operator v„„^ . The second is the scattering 
potential U(r) which can produce interband scattering. 
We quantify the role of interband coupling by separat- 
ing both velocity operators in Eq. <j2j) into intra-band 
and inter-band contributions to express the conductiv- 
ity as the sum of intra-band (oc v±±v±±), inter-band 
(oc v± T v T ±), and interference (oc v±±v± T ) terms. We 



find that the interference terms average to negligible val- 
ues. As illustrated in Fig. [IJinset), the intra-band con- 
tribution dominates in the higher carrier density Boltz- 
mann transport regime, as expected. However, it does 
not completely vanish at zero density as long as the scat- 
tering potential is in play. The inter-band contribution, 
in contrast, increases substantially near the neutrality 
point. Fig. [TJinset) shows that c m i n is due substantially, 
and possibly dominantly, to the non-classical interband 
coherent contribution. 

In an attempt to isolate the source of the peculiar con- 
ductivity behavior we have in Fig. @ compared the nu- 
merical conductivities of our bilayer model with those of a 
decoupled band model in which H — > h 2 (k 2 + ky)<r z /2m. 
The two models have the same density-of-states, but the 
decoupled band model has no interband velocity-operator 
matrix elements, and the scattering potential U (r) is 
not able to couple the bands. The golden-rule relax- 
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ation times of the models are identical when we also let 
u — > uq/i/2 to compensate for the suppression of right- 
angle scattering in the bilayer case. Fig. shows that 
Cmin ~> in the decoupled band model. Deviations from 
the Drude formula at low carrier concentrations in Fig.[2ji 
have a negative sign and are consistent with Anderson in- 
sulator behavior. In Fig. [2] we also see enhanced conduc- 
tivity compared to the Boltzmann model at larger values 
of uq at high carrier densities, which we attribute sim- 
ply to an overestimate of scattering rates by the golden- 
rule expression. The small negative deviation from the 
Boltzmann model at small uq may partially reflect weak 
localizatio n 27 ^ 8 . 

In the zero-temperature limit of the generalized Boltz- 
mann theory, the integrals over wavevector in the ex- 
pression for the interband-coherence conductivity can be 
evaluated to obtain 



Act = — 

2h 



7T 

tan 

2 



(5) 



and the total Boltzmann conductivity will be ctb = 
ctd+Act. It follows that ctb never falls down below 

^min — 

7re 2 /Ah for any choice of parameter values. This value 
agrees with Ref*2£, where a related modified Boltzmann 
approach is combined with a four-band effective Hamil- 
tonian for the carriers, as well as with recent theoreti- 
cal predictions^ 6 , using other closely related approaches. 
Our CT m i n differs from the one obtained for ballistic bi- 
layer graphen o 24 ' 25 i 37 , where the a m \ n is attributed to 
evanescent modes penetrating the sample from contacts. 
We emphasize that Eq. ([5]) should be only seen as the 
rough analytical approximation for our numerical results. 
Eq. ([5]) together with the Drude term fits the numerical 
conductivity curves quite well, but it docs not mean that 
the conductivity minimum is exactly ire 2 /Ah. However, 
the similarity of CT m i n values obtained with different ap- 
proximate approaches might suggest a common underly- 
ing origin in a relationship to the spectral flows associ- 
ated with the topological propertie d 26 ! 38 ' 39 of graphene 
single-layer and bilayer bands. 



To conclude, in our approximate theory the minimum 
conductivity is mainly due to a electric field driven co- 
herence between the conduction and valence bands. Mo- 
mentum space drift due to the electric field does not re- 
populate momenta in a full valence band, as maintained 
in text-book transport theory, but it does drive the sys- 
tem from equilibrium in that it alters the relationship 
between momentum and sublattice pseudospin. There 
is still exactly one electron at each momentum, but the 
momentum states no longer come with definite helicity, 
i.e. are no longer the equilibrium valence band wavefunc- 
tions. As consequence, all valence electrons contribute 
to the conductivity, although the contribution from large 
momenta (O^r > 1) gets suppressed by the larger spin 
precession. Our numerical calculation provide at least 
partial support for this picture. Quantitative discrepan- 
cies might come from not accounting for weak localization 
effects and the influence of disorder on the equilibrium 
state. 

IV. SUMMARY 



We have used numerical exact diagonalization calcula- 
tions to demonstrate i) that the conductivity of bilayer 
graphene in the limit of zero carrier-density <7 m i n ~ e 2 //i, 
ii) that inter-band coherence response plays a key role in 
this property, and iii) that the formation of electron-hole 
puddles due to strong but smooth potential variations 
is not a necessary condition for the minimum conduc- 
tivity phenomena. We believe that our model is rele- 
vant to suspended graphene samples in which charged 
impurities are removed by annealing. When spin and 
valley degeneracy is accounted for we estimate numeri- 
cally cr^jj ~ 8.2 kfl which appears to be consistent with 
current measurements^ -. 
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